Decay constants of P-wave heavy-light mesons from unquenched lattice QCD 
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We review some decays that require knowledge of the decay constants of 0^ heavy-light mesons. 
We compute the decay constants of P-wave heavy-light mesons from unquenched lattice QCD, with 
two degenerate flavours of sea quarks, at a single lattice spacing. The lightest sea quark mass used 
in the calculation is a third of the strange quark mass. For the charm- strange meson we obtain the 
decay constant: fo ^ ~ 340(110) MeV using our normalisation conventions. We obtain the fp^^^"^ 
(static-strange P-wave) decay constant as 302(39) MeV. 

PACS numbers: 11.15.Ha , 12.38.Gc, 14.40.Cs 



I. INTRODUCTION 

The decay constants of heavy-light P-wave mesons 
have a number of important uses in phenomenology [i|,|2|- 
In this paper we use unquenched lattice QCD to com- 
pute the decay constants of 0+ heavy-light mesons, us- 
ing heavy quarks with masses around that of the charm 
quark, and with heavy quarks in the static limit. 

Due to heavy-quark symmetry, there are four P-wave 
[L = 1) heavy-light excited mesons, noted for short D** , 
organised in two doublets, one doublet carrying a total 
angular momentum j = 1/2 and the other j — 3/2 for 
the light quarks. In the charm sector, we use the notation 
that the scalar J^ — Q^{j — 1/2) meson is Dq+ when 
the light quark is an up or a down quark and D^ o+ when 
it is a strange quark. In the static limit of HQET, the 
members of each of the two doublets are degenerate in 
mass. 

The low mass of the Z3s(2317) meson, recently discov- 
ered by BaBar :3] and confirmed by other experiments, 
relative to the quark model predictions of Godfrey and 
Isgur [J, I3 was originally a puzzle. The quantum num- 
bers of Ds{2'in) are consistent with J^ — 0+, but this 
needs confirmation id| . There are some speculations that 
this state may be a molecule ,7|. The theoretical studies 
to understand the DJ22,n) have recently been reviewed 
by Colangelo et al. |3l and by Swanson |^. The lattice 
QCD results for the mass of the Dg 0+ state were tantalis- 
ing close (within large err ors) to the experimental mass 
of the i:»s(2317) state jHllll. However, the lattice cal- 
culation of the Dg 0+ meson in [131 did not agree with the 
mass of the Ds{2in). As discussed by UKQCD T?] a 
hadronic state will couple to many different interpolating 
operators made out of quarks, anti-quarks and glue, with 
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the same quantum numbers as the hadron. In order to 
say whether a hadronic state is more like a molecule or 
qq state, the lattice calculation needs to also determine 
the amplitude for a hadronic state to be in a specific con- 
figuration of quarks and anti-quarks. One possible way 
to determine the quark distribution of a hadronic state 
is to look at observables like form factors, the most basic 
being the decay constant. For example, in a simple pic- 
ture of a molecular state, the decay constant should be 
suppressed relative to that of a bound state of quark and 
antiquark. An analysis of non-leptonic decays of the B 
meson using the factorisation hypothesis suggested that 
the decay constant of the Z^s (2317) could be significantly 
smaller than that of the pseudo-scalar heavy-light me- 
son [ijj ll^ ■ This motivates a lattice calculation of the 
decay constant of the 0+ heavy-light meson. 

The B -^ D**TT decays are also relevant to clarify the 
so called '1/2 vs. 3/2 puzzle' El IH El Indeed, the 
heavy to heavy B — > D** transitions, with j = 1/2, 3/2, 
are parametrised by the generalised Isgur- Wise functions 
Tj{w), where w = vb ■ vd"- Several theoretical con- 
siderations using independent methods (sum rules de- 
rived from QCD [l9l |. covariant quark model [23, lattice 
QCD |2lL 1221 ■ e xperimental data combined with naive 
factorisation E3)i suggest that the j — 3/2 states dom- 
inate, with respect to the j — 1/2 states in _B ^ D** 
transitions, i.e. T3/2(1) > ti/2(1). Nevertheless this 
prediction is in contradiction with some of the currently 
available experimental data. The solution for this puzzle 
still requires new inputs from both theory and experi- 
ment. From the theoretical side, the determination of 
the decay constant, fu + of the Dq+ meson is needed 
to evaluate the decay width of two (out of the three) 
classes of B -^ D**7r decays, when naive factorisation is 
assumed |17l] . 

A recent summary [l3| of results for the decay con- 
stant, showed that estimates for fu ^ lay between 122 
and 417 MeV. This summary also reported that the de- 
cay constant of 0+ heavy-light mesons in the static limit 



was higher than at the charm mass. This required con- 
firmation, because the heavy quark mass dependence 
was obscured by model dependence. Some factorisation 
schemes zJ| for non-leptonic decays of B mesons use the 
static hniit as the leading term, hence computing the 
decay constant of the P-wave static-light meson is also 
important for this reason. The evaluation of the l/mg 
corrections to the heavy quark limit of the various quan- 
tities describing B -^ D** transitions is crucial since 
simple estimates tend to suggest the presence of large 
corrections [P^ . 

There is a huge effort in lattice gauge theory that aims 
to compute the decay constant of the heavy- light pseudo- 
scalar mesons, because they are crucial for determining 
CKM matrix elements |2J, |23 . The decay constant of 
the heavy-light vector meson has been computed in |2q| 
in order to test the HQET scaling relations. 

There has been no published work (apart from a pri- 
vate communication reported in [l7J) on using lattice 
QCD to compute the decay constants of P-wave heavy- 
light mesons. The study of the wave-function of P-wave 
static-light mesons 27], using quenched lattice QCD, did 
not extract a decay constant. UKQCD has studied the 
charge and matter radial distributions of P-wave static- 
light mesons [2ll29l|. 

In our calculations, we used the non-perturbatively im- 
proved clover action with clover coefficient of 2.0171 for 
the light quarks. The Wilson gauge action was used. The 
lattice volume was 16^ 32 and (3 = 5.2. There were two 
flavours of sea quarks. The k values are tabulated with 
the results. The sea quark masses span from the strange 
quark to one third of the strange quark mass. The details 
of the light hadron spectroscopy can be found in 30, 31| ■ 

The plan of the paper is that in section ^ we first dis- 
cuss our results using clover fermions with quark masses 
around the charm value. Next in section ITTTl we present 
data using static heavy quarks. In the final section Hvl 
we compare our results to other determinations and make 
some remarks about the heavy quark mass dependence 
of the decay constants. 



II. DECAY CONSTANTS AROUND THE 
CHARM MASS 

We have used the unqucnched data from [l^, |32| to 
study the decay constant of the 0"*" meson around the 
charm mass with clover quarks. 

The data set used the sea k = 0.1350 and three va- 
lence light quark masses using K-values of 0.1340, 0.1345 
and 0.1350. The k values for the heavy quarks were: 
0.113, 0.119 and 0.125. Usingro = 0.5 fm, the inverse 
lattice spacing is 1.88 GeV |30|- The sea quark mass 
was kept fixed at the mass of the strange quark, hence 
this is a partially quenched calculation. The data sample 
included 394 gauge configurations separated by 20 trajec- 
tories. Adjacent data were binned together. A smearing 
matrix of order two, using local and fuzzed basis states, 



was fitted to a factorising fit model |30j . 

The decay constant of the Dq+ meson can be defined 
by equation. 
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where V^'^ is the vector current with two non-degenerate 
quark flavours. This uses the V — A structure of the 
charged weak quark current and parity conservation. We 
used a unit gamma matrix between the heavy and light 
quark fields as the interpolating operator for the Dq+ 
state. Unfortunately, we did not measure the lattice cor- 
relators with 7o at the sink and 1 at the source that would 
be required to use equation ^ directly. In principle due 
to the number of gamma matrices, there are 256 different 
local meson correlators. Many lattice QCD codes do not 
compute every possible correlator, because some corre- 
lators are zero by symmetries. The correlator with a 70 
at the sink and a unit matrix at the source is zero for 
degenerate quarks because of charge conjugation. 

Colangelo et al. [l| define the decay constant {go+) 
using equation 13 
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where rric is the mass of the charm quark. Equation [21 
can be related to equation Q] since the divergence of the 
vector current for non-degenerate quarks is proportional 
to the scalar density. The derivation of equation [21 from 
equation ^assumes that the mass of the light quark can 
be neglected relative to that of the charm quark. 

From a lattice QCD perspective, it is more natural to 
define the decay constant using equation [3 



(0|cg|Do+)=Mn ,/o+ 



(3) 



The determination of the mass of the charm quark can 
have quite large systematic errors at a fixed lattice spac- 
ing |32| , so we prefer the definition of the decay constant 
with no explicit factor of the charm mass. To convert 
from our normalisation of the decay constant /o+ to that 
of gQ+ , we note that for this data set UKQCD obtained 

mdmcY^ = 1.247(3)lf GeV [al, using a quark mass 
defined using the Fermilab heavy quark formalism ^53]. 

The matrix element is related to the coupling ampli- 
tude Zi in the fit to the 0+ to 0+ meson correlator. 
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A simple linear fit model was used to extrapolate and 
interpolate in the quark masses. We first interpolated 
to the strange quark mass 3^, or extrapolated to the 
light quark mass (at Kcrit)- The decay constants were 
then interpolated to the value of the mass of the charm 
quark 32] . 

In figure^we report an effective mass plot for the local- 
local and fuzzed-local correlators used in the analysis. In 
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FIG. 1: Effective mass plot and fitted model for the correla- 
tors with heavy k — 0.113 and light k = 0.1350. 



FIG. 3: Mass dependence of the P-wave decay constant de- 
fined in equation 13 at the strange quark mass, as a function 
of the bare heavy quark mass in physical units. The vertical 
line is the estimate for the charm quark mass. 



O data 

□ at strange 



FIG. 2: Mass dependence of the P-wave decay constant de- 
fined in equation 13 at the fixed heavy quark mass with k — 
0.113, as a function of the bare light quark mass in physical 
units. The estimate of the decay constant at the mass of the 
strange quark is also shown. 



figure 121 we plot the bare decay constant as a function of 
the bare light quark mass in physical units with mass of 
the heavy quark fixed. Similarly, figure 13 shows the bare 
lattice decay constant as a function of the heavy quark 
mass, with the light quark mass fixed at the mass of the 
strange quark. The dependence of the decay constant on 
the mass of the quarks is very mild. We discuss some of 
the issues in the chiral extrapolations in section UTTl 



To convert the lattice number into the MS scheme we 
use tadpole improved perturbation theory to one loop 
order (the required expressions are listed in J35ll36l |). 



Zs = uo{l 
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where uq is the fourth root of the plaquette. The mass 
dependent improvement factor (1 + bg ™' 2™^ ) was mul- 
tiplied into the scalar current. We used the one loop 
expression for bs |37j . 



6s = (l + aj.3722) 



(6) 



and the same prescription for the coupling as used in 
UKQCD's estimate of the mass of the charm quark |32 |. 



We use /I 



to determine Zs- 



In lattice units we obtain fo ^/Zs = 0.280(70) and 
fD^^+/Zs = 0.270(90). Using ro = 0.5 fm, and Zs from 
equation we obtain fo + = 360(90) MeV and fo + 
— 340(110) MeV. In principle, we should estimate the 
effect of the lattice spacing dependence and higher order 
terms in the perturbative expansion, as we did when we 
estimated the mass of the bottom quark J23|. However, 
because the statistical errors are of the order of 30%, the 
statistical errors will clearly dominate. For example, a 
typical estimate of the systematic error due to the differ- 
ent ways of determining the lattice spacing in a lattice 
calculation with these parameters would be to vary the 
value of ro between 0.5 fm and 0.55 fm. This systematic 
error is a 10 % effect. 

We now discuss the central value of ro used. When 
To, a number derived from the heavy quark potential, 
was introduced by Sommer [39j, he quoted a value of 
To about 0.49 fm. The value of ro can also be deter- 
mined from lattice calculations. As reviewed in [l(]| . un- 
quenched calculations with similar parameters to those 
used here compute a value of ro between 0.5 and 0.55 



fm. The issue is not so clear because the HPQCD col- 
laboration, using gauge configurations from the MILC 
collaboration, obtain a value of rp — 0.469(7) fm. As 
UKQCD has argued [Sg, for calculations with these pa- 
rameters it better to use a value of tq closer to 0.5 fm 
because this is the value determined from similar calcu- 
lations. The MILC/HPQCD value for tq is with in 10% 
of 0.5 fm. 



III. DECAY CONSTANTS IN THE STATIC 
LIMIT 

This part of the work is a continuation of UKQCD 's 
study of static- light mesons. The static-light meson spec- 
trum was reported in [l^, and then used to extract the 
mass of the bottom quark |23 . Here we use the extended 
data set that was used to look for chiral logs in the fs 
decay constant |4nl |. 

In addition to using the standard static action of 
Eichten and Hill, we also report data using one of the 
static actions developed by the ALPHA collaboration |4y 
with improved signal to noise ratio. The exact static ac- 
tions are described in |40| . 

Here we focus on the meson doublet with total light 
quark momentum of j = 1/2, we call this the P_ static- 
light meson [l3|. As corrections to the static limit are 
included the P_ state will split into a 0+ state and a 1+ 
state. In this paper, the decay constant of the P_ static- 
light meson is called /|*°*"=. The /|*°*"= decay constant 
is defined by the matrix element in equation 
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where V^t is the vector current [l|, |2|. The simple spin 
structure of the static quark means the same decay con- 
stant is obtained if in equation|7|the vector current (with 
jjL = 0) \s replaced with a scalar current. The /|*''**'= ma- 
trix element is extracted from the amplitudes in the two 
point correlator 

C{t) = ^(0|Fo(x,t)*tj(x,0)|0) (8) 

X 

^ ^£*«tic^^^ exp{~a£t) , (9) 

where ^^b is the interpolating operator for the P- static- 
light meson and ground state dominance is shown in 
equation El In practise we use a basis of smearing func- 
tions to do variational smearing 13]. The 2'£*''**'^ ampli- 
tude is related to the /|*°*"= decay constant 



Mb. 



(10) 



where Zp^"'*'"'^ is a perturbative matching factor that we 
discuss below. This is an equivalent definition to the 
one used for the pseudo-scalar heavy-light decay constant 

/ J'static\ 



It is traditional to match the results of a lattice static- 
light calculation to continuum QCD via two steps p3 |. 
QCD is matched to the continuum static theory. The lat- 
tice static theory is then matched to the static continuum 
theory. 

The matching of the static continuum theory to the 
static-lattice theory has been done by Eichten and 
Hill J43 and by Borrelli and Pittori 44] for the clover 
action. This matching was done for a heavy-light current 
with arbitrary gamma matrix (P). For the definition of 
the decay constant in equation we do the perturbative 
matching assuming P = 70. The matching between the 
continuum static theory and the static lattice theory is 
via Z(P). 



Z{T) = 1 
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where 



Ar = di + {d2 - d^)G ^ 
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(11) 
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where the values of the constants are di — 5.46, d2 = 
-7.22, / = 13.35, e = 4.53 0, d^ = -4.04, and /^ = 
-3.63 |4^. G = 1 when P = 70. 

The matching of the continuum static theory to con- 
tinuum QCD was done using [4J] 



= 1 
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where mb is the mass of the bottom quark. The equiv- 
alent expression for equation I13L for arbitrary P matrix, 
is in [ijl ■ The matching factor Zf"-^^'^ in equation is the 
product of Z{T = 70) with Zg^tat to leading order in the 
square of the coupling g^. 

The improvement coefficients have not been calculated 
for the P-wave decay constant and we will therefore not 
include them in our analysis. We used a simple boosted 
couphng (g^/(uo)^ and g^ = |) to compute the pertur- 
bative matching factors. 

We use the mass (5279-1-400 MeV) for the mass factor 
in eauation llOl The experimental spectrum of the P-wave 
B mesons is not very well determined at the moment |6|, 
I3. Lattice calculations predict that the lowest P-wave 
meson will be roughly 400 MeV above the S-wave states. 
So we add 400 MeV to the mass of the B+ state 0. 
This corresponds to a 4% effect on the decay constant. 
We prefer not to use the HQET scaling law of the decay 
constant to quote numbers for the decay constant of the 
0+ charm-light meson, as is done by Jugeau et al. [l7| . 
because it is known that 1/M corrections to the static 
limit of the decay constant of the pseudo-scalar heavy- 
light meson are large. 

The unrenormalised data for the amplitudes are in ta- 



ble |III| The final column shows the value for J|*°**'= in 
the static limit in physical units. The value of rg = 0.5 
fm was used to convert the data into physical units. 



Name 


formalism Ksea ^J**"**"^ 


/!,*«*- MeV 


DF3 


Eichten-Hill 0.1350 0.240^27 


302(39) 


DF3 


Fuzzed ALPHA 0.1350 0.199+20 




DF4 


Eichten-Hill 0.1355 0.234ti^ 


322(59) 


DF4 


Fuzzed ALPHA 0.1355 0.174+^° 




DF6 


Eichten-Hill 0.1358 0.182^36 


272(78) 


DF6 


Fuzzed ALPHA 0.1358 O.I31I35 





TABLE I: Amplitudes and decay constants for the P_ static- 
light mesons. Further information about the lattice parame- 
ters is in Hol. 



As the renormalisation factor for the ALPHA static 
action has not yet been determined, we can't use that 
data to quote a physical number. What is disappoint- 
ing is that the ALPHA static action does not produce 
any reduction in the statistical errors over the standard 
Eichten-Hill static action in our case. For the static-light 
pseudo-scalar meson decay constant the statistical errors 
were significantly smaller for the ALPHA static action 
than the Eichten-Hill action [40| . 

In [lj| it was shown that the DF3 data set corre- 
sponded to sea quarks with mass around the strange 
quark mass. Hence we quote the /|*''*"= decay constant 
as 302(39) MeV. 

The value of the P-wave 0+ decay constant at the 
strange quark mass is of interest to compare to model 
calculations or for future decays of the Bs meson mea- 
sured at a hadronic experiment such as LHCb, CDF, or 
DO. The decay constant required for the factorisation 
analysis of the heavy to heavy B -^ D**tt is the /|*''*"= 
decay constant in the chiral limit |17| . 

To study the light P-wave 0+ decay constant the chi- 
ral extrapolation must be discussed. The importance of 
chiral logs in the mass extrapolation of the Jb decay con- 
stant has only recently been observed as the dominant 
systematic error in the determination of fs/ fs^ |45ll4g|. 

The equivalent chiral perturbation theory calculation 
for the P-wave decay constant (to our knowledge) has 
not yet been done. There are calculations of the masses 
of heavy-light mesons to one loop in heavy-light chiral 
perturbation theory [JTI . l4q . UKQCD has previously 
computed the relevant hadronic coupling [43 • This result 
is compared against experiment in |47l |. 

The data for the P-wave decay constant in table lHll are 
not really precise enough to determine the quark mass 
dependence. A simple linear fit against the square of the 
pion mass of the data in table UTTI gives /|*°**'= = 294(88) 
MeV in the chiral limit of zero mass light quarks. 



IV. CONCLUSION 

We have computed the decay constant of 0+ heavy- 
light mesons using an unquenched lattice QCD calcula- 
tion at a single lattice spacing. For the static-light decay 



Method 


go+ MeV 


This work 


200 ± 50 


QCD sum rules [Ij 


170 ± 20 


Lattice QCD [Hi 


122 ± 43 


B -^ D**-K fl7] 


206 ± 120 



TABLE H: Comparison of decay constants of charm-light 
Dq+ meson. 



constant, we used sea quarks as low at a third of the 
strange quark mass. For the calculation of the 0+ decay 
constant, we did a partially quenched analysis with the 
sea quark fixed at the strange quark mass. 

We obtain the /|>f *'== decay constant as 302(39) MeV. 
Given the qualifications mentioned in the previous sec- 
tion, we obtain /|,*"**'= = 294(88) MeV. The data for 
jpatic fj-Qj-fj various models are summarised in jTSL Two 
QCD sum rule results are /l,*"'*'^ =304 ± 400 and 
377 ± 53 MeV j^ . Using clover quarks for the charm 
quark, we obtain fo + ^ 360(90) MeV and fo + = 
340(110) MeV. The magnitude of these decay constants 
can be compared against the value of the fo, and fo 
pseudo-scalar decay constants. From unquenched lat- 
tice QCD, the Fermilab, MILC and HPQCD collabo- 
rations 51] obtained: fo, = 249 ± 3 ± 16 MeV and 
/£,+ = 201 ± 3 ± 17 MeV. The CLEO-c collabora- 
tion have recently reported |52l | the experimental result 
/d+ = 223(17) (3) MeV. 

Although the systematic errors are different in the two 
lattice calculations, we can make the qualitative state- 
ment that the P-wave and S-wave decay constants of 
charm-light mesons are similar in magnitude. This cal- 
culation does not support a suppression of the decay con- 
stant of the P-wave heavy-light meson relative to the de- 
cay constant of the S-wave heavy-light meson. 

The only decay constant of a light P-wave meson that 
is sometimes calculated from lattice QCD is that of the 
ai meson. The decay constant of the ai meson, denoted 
as fai , is usually computed from a definition that has 



dimension MeV . If we use the results from Wingate 
et al. [Ill, then the value of ^ is 240 MeV. This is 
larger than the value of the pion decay constant 131 MeV, 
showing a similar trend to the heavy-light case. 

Kurth and Sommer [5j have noted that there are po- 
tential problems with extrapolating in the heavy-quark 
mass unless the continuum limit has been taken. Also 
see Kronfeld 55] for a discussion of the problems com- 
bining static and propagating heavy-light data. Given 
these theoretical concerns, the data is consistent with the 
view that the decay constant of the static-light meson is 
larger than that of the 0+ charm-light meson, if we use 
definition of the decay constant in equation ^ (our value 
of 5o is in table 2) , but the errors need to be reduced for 
a definitive statement. 

Jugeau et al. [i3|. have collected together a number 
of different calculations of the 0+ decay constant of a 



charm-light meson. To compare our results to other cal- 
culations we use the normalisation in equation ^ Hence 
we multiply our result by TOc/-^Ds(2317) — 1.27/2.317 and 
compare to other calculations in table 2. There is rea- 
sonable agreement between different determinations. To 
reduce the errors on the lattice results requires a calcula- 
tion of similar effort to that done to compute the decay 
constants of the heavy-light pseudo-scalar mesons |51(| . 



supported by EPSRC grant GR/K41663, PPARC grants 
GR/L22744 and PPA/G/S/ 1998/00777. This work has 
been supported in part by the EU Integrated Infrastruc- 
ture Initiative Hadron Physics (I3HP) under contract 
RII3-CT-2004-506078. We are grateful to the ULgrid 
project of the University of Liverpool for computer time. 
One of the authors (CM) wishes to thank PPARC for the 
award of a Senior Fellowship. 



V. ACKNOWLEDGEMENTS 

We thank Alex Dougall for help. The gauge configura- 
tions were generated on the Cray T3E system at EPCC 



[1] P. Colangelo, G. Nardulli, A. A. Ovchinnikov, and 

N. Paver, Phys. Lett. B269, 201 (1991), 
[2] P. Colangelo, G. Nardulli, and N. Paver, Phys. Lett. 

B293, 207 (1992), 
[3] BABA R, B. Au bert et ai, Phys. Rev. Lett. 90, 242001 

(2003), |hep-ex/ 0304021 
[4] S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985), 
[5] S. Godfrey and R. Kokoski, Phys. Rev. D43, 1679 

(1991), 
[6] Particle Data Group, S. Eidelman et ai, Phys. Lett. 

B592, 1 (2004), 
[7] T. Barnes, F. E. Close, and H. J. Lipkin, Phys. Rev. 

D68, 054006 (2003), hep-ph/0305025 
[8] P. Colangelo, F. De Fazio, and R. Ferrandes, Mod. Phys. 

Lett. A19, 2083 (2004), hep-ph/0407137 
[9] E. S. Swanson, (2006), hep-ph/0601110 
[10] UKQCD, A. Dougall, R. D. Kenway, C. M. May- 
nard, and C. McNeile, Phys. Lett. B569, 41 (2003), 
hep-lat/0307001j 
[11] RBC, S. Ohta, H. Lin, and N . Yamada, PoS LAT2005, 

096 (2005), hep-lat/0510071 
[12] G. S. Bali, Phys. Rev. 

hep-ph/0305209 
[13] UKQCD, A. M. Green, J. 
C. Michael, and G. Thompson, 
[14] A. Datta and P. J. O'donnell, 

(2003), hep-ph/0307106 
[15] D. S. Hwang and D.-W. Kim, 

(2005), hep-ph/0410301 
[16] N. Uraltsev, (2004), hep-ph/0409125| 
[17] F. Jugeau, A. Le Yaouanc, L. Oliver, and J. C. Raynal, 

Phys. Rev. D72, 094010 (2005), hep-ph/0504206 
[18] L L Bigi et al, (2005), hep-ph/0512270 
[19] N. Uraltsev, Phys. Lett. B501, 86 (2001), 

hep-ph/0011124 
[20] A. Le Yaouanc, L. Oliver, O. Pene, and J. C. Raynal, 

Phys. Lett. B365, 319 (1996), hep-ph/9507342 
[21] D. Becirevic et ai, Phys. Lett. B609, 298 (2005), 

hep-lat/0406031 
[22] B. Blossier, A. Le Yaouanc, V. Morenas, and O. Pene, 

Phys. Lett. B632, 319 (2006), hep-lat/0507024 
[23] M. J. Dugan and B. Grinstein, Phys. Lett. B255, 583 

(1991), 
[24] V. Lubicz, Nucl. Phys. Proc. Suppl. 140, 48 (2005), 



D68, 071501 (2003), 

Koponen, C. McNeile, 
(2003), hep-lat/0312007 
Phys. Lett. B572, 164 

Phys. Lett. B606, 116 



[25; 
[26; 

[27 

[28; 

[29 
[30 
[31 

[32; 

[33 

[34; 
[35; 
[36 
[37 
[38; 
[39 

[4o; 

[41 

[42' 
[43^ 
[44' 
[45| 

[46; 

[47 



|hep-lat/041005T] 

M. Okamoto, PoS LAT2005, 013 (2005), 

hep-lat/0510113 

C. Bernard et ai, Phys. Rev. D65, 014510 (2002), 
hep-lat/0109015 

A. Duncan, E. Eichten, and H. Thacker, Phys. Lett. 

B303, 109 (1993), 

UKQCD, A. M. Green, J. Koponen, and C. Michael, AIP 

Conf. Proc. 756, 369 (2005), hep-lat/0412002 

UKQCD, A. M. Green et ai, PoS LAT2005, 205 (2005), 

hep-lat/0509161 

UKQCD, C. R. AUton et ai, Phys. Rev. D65, 054502 

(2002), hep-lat/0107021 

UKQCD, C. R. AUton et ai, Phys. Rev. D70, 014501 

(2004), hep-lat/0403007 

A. Dougall, C. M. Maynard, and C. McNeile, JHEP 01, 

171 (2006), hep-lat/0508033 

A. X. El-Khadra, A. S. Kronfeld, a nd P. B. M ackenzie, 

Phys. Rev. D55, 3933 (1997), hep-lat/9604004' 

UKQCD, D. J. Hepburn, Nucl. Phys. Proc. Suppl. 119, 

338 (2003), hep-lat/0209036 

M. Gockeler et al, Nucl. Phys. Proc. Suppl. 53, 896 

(1997), hep-lat/9608033 

T. Bhattacharya, R. Gupta, W.-J. Lee, and S. R. Sharpe, 

Phys. Rev. D63, 074505 (2001), hep-lat/0009038 

S. Sint and P. Weisz, Nucl. Phys. Proc. Suppl. 63, 856 

(1998), hep-lat/9709096 

UKQCD, C. McNeile, C. Michael, and G. Thompson, 

Phys. Lett. B600, 77 (2004), hep-lat/0408025 

R. Sommer, Nucl. Phys. B411, 839 (1994), 

hep-lat/9310022 

UKQCD, C. McNeile and C. Michael, JHEP 01, Oil 

(2005), hep-lat/0411014 

ALPHA, M. Delia Morte et ai, Phys. Lett. B581, 93 

(2004), hep-lat/0307021 

E. Eichten and B. Hill, Phys. Lett. B234, 511 (1990), 

E. Eichten and B. Hill, Phys. Lett. B240, 193 (1990), 

A. Borrelli and C. Pittori, Nucl. Phys. B385, 502 (1992), 

JLQCD, S. Aoki et ai, Phys. Rev. Lett. 91, 212001 

(2003), hep-ph/0307039, 

A. S. Kronfeld and S. M. Ryan, Phys. Lett. B543, 59 

(2002), hep-ph/0206058 

D. Becirevic, S. Fajfer, and S. Prelovsek, Phys. Lett. 
B599, 55 (2004), hep-ph/0406296 



[48] T. Mehen and R. P. Springer, Phys. Rev. D72, 034006 

(2005), hep-ph/0503134 
[49] UKQCD, C. McNeile, C. Mich ael, and G. Thom pson, 

Phys. Rev. D70, 054501 (2004) , |hep-lat/0404010' 
[50] S.-L. Zhu and Y.-B. Dai, Mod.Thys. Lett. A14, 2367 

(1999), hep-ph/9811449 
[51] C. A ubin et al, Phys. Rev. Lett. 95, 122002 (2005), 
|hep4a t/0506030i 



[52] CLEO, M. Artuso et al, Phys. Rev. Lett. 95, 251801 

(2005), hep-ex/0508057, 
[53] M. Wingate, T. DeGrand, S. Collins, and U. M. Heller, 

Phys. Rev. Lett. 74, 4596 (1995), hep-ph/9502274 
[54] ALPHA, M. Kurth and R. Sommer, Nucl. Phys. B623, 

271 (2002), hep-lat/0108018 
[55] A. S. Kronfeld, (2002), hep-lat/020502l| 



